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Composite operators in lattice QCD: nonperturbative renormalization 
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We investigate the nonperturbative renormalization of composite operators in lattice QCD restricting ourselves 
to operators that are bilinear in the quark fields. These include operators which are relevant to the calculation 
of moments of hadronic structure functions. The computations are based on Monte Carlo simulations using 
quenched Wilson fermions. 



1. INTRODUCTION 

Lattice investigations of hadronic structure re- 
quire the calculation of hadronic matrix elements 
of composite operators. In general, one has to 
convert the bare lattice operators into renor- 
malized continuum operators by multiplication 
with the appropriate renormalization coefficient 
Z. The accuracy of perturbative calculations of 
these Z factors remains uncertain, even if tadpole 
improvement |l| is used, as more than a one-loop 
computation is rarely available. Therefore a non- 
perturbative calculation by Monte Carlo simula- 
tions appears to be an attractive alternative . 

We have performed such a computation for 
operators that determine moments of hadronic 
structure functions. Computing the Z factors for 
a rather large range of renormalization scales we 
try to find out at which scales (if at all) pertur- 
bative behavior sets in such that the multiplica- 
tion with perturbative Wilson coefficients makes 
sense. For further details and results see Ref. || . 

We use standard Wilson fermions in the 
quenched approximation. At = 6.0 ((3 — 6.2) 



we work on a 16 3 x 32 (24 3 x 48) lattice. For both 
(3's we have studied three values of the hopping 
parameter k. 

In this talk, we shall concentrate on one par- 
ticular example, namely the operator 



= q 
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(0 V2 b in the notation of Ref. ||) whose hadronic 
matrix elements are proportional to the momen- 
tum fraction (x) carried by the quarks. Work- 
ing with Wilson fermions it is straightforward to 
write down a lattice version of the above opera- 
tor. One simply replaces the continuum covariant 
derivative by its lattice analogue. 

2. THE METHOD 

We follow the procedure proposed by Martinelli 
et al. |2j. It applies the definitions used in (con- 
tinuum) perturbation theory to vertex functions 
calculated by Monte Carlo simulations on the lat- 
tice. To be more specific, we calculate the quark- 
quark Green function (a matrix in color and Dirac 
space) with one insertion of the operator O at 
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momentum zero in the Landau gauge. Amputat- 
ing the external quark legs we obtain the vertex 
function T(p) as a function of the quark momen- 
tum p. Defining the renormalized vertex function 
by = Z~ 1 ZT(p) we fix the renormalization 

constant Z by imposing the condition 



£tr (rR^r^tp)- 1 ) = 1 



(2) 



at p 2 = /i 2 , where /i is the renormalization scale. 



So we calculate Z from 



Z x Z^ix (r(p)r Bo rn(p)~ 
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(3) 



with p 2 = p? . Here rBorn(p) is the Born term in 
the vertex function of O computed on the lattice, 
and Z q denotes the quark field renormalization 
constant. The latter is calculated from the quark 
propagator S[p): 



Z q (p) 



(4) 



tr (-iJ2\l* s ' m ( a P\) aS Hp)) 
12 £ A sin 2 (ap A ) 

again at p 2 = fi 2 . 

3. PERTURBATIVE INPUT 



Eq. (||) defines a renormalization scheme of the 
momentum subtraction type, which we call MOM 
scheme. In order to convert our results into the 
more popular MS scheme we have to perform a 
finite renormalization. The corresponding renor- 
malization constant Z^q M is computed in contin- 
uum perturbation theory using dimensional reg- 
ularization. 

Neglecting quark masses we obtain in the Lan- 
dau gauge 



7 MS 
J MOM 
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pi))' 



3p 2 (pl + ±(pl+pi+pt)) 



(5) 



0{g A ). 



where Cf — 4/3 for the gauge group SU(3). Be- 
cause of the noncovariance of the condition (0), 
i^MOM depends on the direction of the momentum 
P- 



For the scale dependence, the renormalization 
group predicts (at fixed bare parameters) 



R(n,Hn) := 



exp 




(6) 



in terms of the running coupling <?(/^ 2 ), the 8- 
function 0(g), and the anomalous dimension j(g). 
An analogous formula describes the 8 dependence 
at fixed renormalized quantities. 

4. RESULTS 

Calculating the necessary two- and three-point 
functions with the help of momentum sources 
P we achieve small statistical errors with a mod- 
erate number of configurations. After extrapolat- 
ing to the chiral limit (linearly in 1/k) we multiply 

by -^mom- 

For the presentation of our results we convert 
lattice units to physical units using a -2 = 3.8 
GeV 2 [p = 6.0) and a~ 2 = 7.0 GeV 2 (8 = 6.2) 
as determined from the string tension p| . The A 
parameter is taken to be A^jg = 230 MeV. 

The scale dependence of our Z's should be de- 
scribed by the renormalization group factor R 
(see @). Hence we divide our numerical re- 
sults by this expression (evaluated in two-loop 
approximation with Hq = 4 GeV 2 ) and define 
Zrgi = Z/R. For Zrgi we hope to obtain a [i 
independent answer, at least in a reasonable win- 
dow of [i values. There /i 2 would be large enough 
to allow for perturbative scaling behavior. On the 
other hand, /i 2 should be small enough to avoid 
strong cut-off effects. 

In Fig. |l| we plot Zrgi versus the renormaliza- 
tion scale pi 2 for our operator O. The errors are 
purely statistical. A "flat" region seems to start 
at [i 2 w 5 GeV 2 . Note that the = 6.0 data 
and the (perturbatively rescaled) 3 = 6.2 data 
agree except for the lowest values of /i 2 . This indi- 
cates that the observed \i dependence is physical. 
For other operators, the results look qualitatively 
similar, although in most cases a window starts 
only at /i 2 w 10 GeV 2 , where it is hard to believe 
that cut-off effects are negligible. 
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Figure 1. Zrgi for the operator O (Eq. (§)). The 
open circles (filled squares) represent the data for 
(3 = 6.2 (/3 = 6.0). The (3 = 6.2 data have been 
rescaled perturbatively to /3 = 6.0. 
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Figure 2. Collection of Z's in the MS scheme 
versus the number of derivatives in the opera- 
tors. The squares represent the nonperturbative 
results, the diamonds (filled circles) correspond to 
the values obtained in lattice perturbation the- 
ory (with tadpole improvement). The f3 — 6.0 
(/3 = 6.2) data are displaced to the left (right). 

In Fig. ^| we compare the nonperturbative val- 
ues for Z in the MS scheme with (tadpole im- 
proved) one-loop lattice perturbation theory. We 
plot the results for all operators with derivatives 
studied in Ref. H at /i 2 = a -2 versus no, the 
number of derivatives. One-loop perturbation 
theory underestimates the increase of Z with nu- 
Tadpole improvement works in the right direction 
for no > 1 without being quantitatively satisfac- 
tory. However, the differences between the three 
methods decrease as j3 (and /i 2 ) becomes larger. 



5. SUMMARY 

In summary, we have performed a compre- 
hensive study of nonperturbative renormalization 
for various operators in the framework of lat- 
tice QCD. Twist-2 operators appearing in unpo- 
larized (polarized) deep-inelastic scattering have 
been studied for all spins < 4 (< 3) j|. The scale 
dependence of Z for the operator considered in 
this talk agrees with perturbative expectations to 
a good approximation already for moderate val- 
ues of /x 2 . However, most of the other operators 
containing covariant derivatives seem to approach 
perturbative scaling only for rather large values 
of the renormalization scale where cut-off effects 
might contaminate the results. The behavior at 
smaller /i 2 could be nonperturbative physics, but 
at present we cannot exclude the possibility that 
it is faked, e.g., by finite-size effects. In any case, 
it might not always be sensible to combine these 
nonperturbative Z's with perturbative Wilson co- 
efficients using a scale of, say, /i 2 = 4 GeV 2 . We 
have therefore started a nonperturbative calcula- 
tion of the Wilson coefficients. First results are 
reported in Ref. 0. 
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